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Antideuterons are among the most promising galactic cosmic ray-related targets for dark matter 
indirect detection. Currently only upper limits exist on the flux, but the development of new experi- 
ments, such as GAPS and AMS-02, provides exciting perspectives for a positive measurement in the 
near future. In this Paper, we present a novel and updated calculation of both the secondary and 
primary d fluxes. We employ a two-zone diffusion model which successfully reproduces cosmic-ray 
nuclear data and the observed antiproton flux. We review the nuclear and astrophysical uncertain- 
ties and provide an up to date secondary {i.e. background) antideuteron flux. The primary {i.e. 
signal) contribution is calculated for generic WIMPs annihilating in the galactic halo: we explicitly 
consider and quantify the various sources of uncertainty in the theoretical evaluations. Propagation 
uncertainties, as is the case of antiprotons, are sizeable. Nevertheless, antideuterons offer an exciting 
target for indirect dark matter detection for low and intermediate mass WIMP dark matter. We 
then show the reaching capabilities of the future experiments for neutralino dark matter in a variety 
of supersymmetric models. 

PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,12.60. Jv,95.30.Cq 



I. INTRODUCTION 

The identification and the understanding of the nature 
of dark matter (DM) is one of the deepest open problems, 
together with the solution to the dark energy mystery, 
in the fundamental physics research. Many experimen- 
tal efforts devoted to the detection of the astronomical 
dark matter in the halo of our and nearby galaxies have 
been carried out in underground laboratories, in large- 
area surface telescopes as well as in space. In the near 
future, the LHC will provide us with invaluable informa- 
tion on particle physics extending beyond the Standard 
Model, thus probing a wide class of theoretical models 
hosting the most viable DM candidates. 

The indirect dark matter detection is based on the 
search for anomalous components due to the annihila- 
tion of DM pairs in the galactic halo, in addition to the 
standard astrophysical production of neutrinos, gamma 
rays and light antimatter in cosmic rays. Data on ncutri- 
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nos, gamma rays, positrons and antiprotons are already 
available at a sensitivity level allowing some inspection 
on possible exotic contributions. In the seminal Paper 
[lj], it was proposed to look for cosmic antideuterons (d) 
as a possible indirect signature for galactic dark mat- 
ter. It was shown that the antideuteron spectra deriv- 
ing from DM annihilation is expected to be much flat- 
ter than the standard astrophysical component at low 
kinetic energies, T-g < 2-3 GeV/n. This argument mo- 
tivated the proposal of a new space-borne experiment 
looking for cosmic antimatter (antiproton and an- 
tideuteron) and having the potential to discriminate be- 
tween standard and exotic components for a wide range 
of DM models. AMS-02 has also interesting capabilities 
of looking for cosmic antideuterons 0]. Antideuterons 
have not been measured so far, and the present experi- 
mental upper limit Q is still far from the expectations 
on the secondary antideuteron flux which are produced 
by spallation of cosmic rays on the interstellar medium 
[7|, l8[ , but in fact perspectives for the near future are very 
encouraging. 

In the present Paper we update and improve our calcu- 
lation of the antideuteron primary flux in a full two-zone 
diffusion model, consistent with a number of indepen- 



2 



dent cosmic ray (CR) measurements, and explicitly esti- 
mate the uncertainties which affect the signal determina- 
tion. In addition, we also provide a new determination 
of the secondary component. In Sect. [Til the framework 
and ingredients for the solution of the two-zone transport 
equation are recalled. In Sect. IIII1 the coalescence model 
for the nuclear fusion process is discussed both for sec- 
ondary and primary d production. Sect. ITVl is dedicated 
to the secondary d flux and to the possible uncertain- 
ties affecting its evaluation. In Sect. [V] the production 
of antideuterons from DM particles is detailed and re- 
sults on the propagated fluxes are presented, together 
with the estimation of the uncertainties due to propa- 
gation, to the dark matter halo profile and to the DM 
annihilation final states. Finally Sect. IVll demonstrates 
that antideuterons are probably one of the most power- 
ful dark matter indirect detection channel, and shows the 
optimistic potentials of next-to-come balloon and space 
based missions. We finally draw our conclusions in Sect. 



II. THE PROPAGATION MODEL 

Cosmic ray fluxes are determined by the transport 
equation as given, e.g., in Berezinskii et al. If steady- 
state is assumed, the transport equation for any nuclear 
species can be rewritten in terms of the cosmic ray dif- 
ferential density dn(r)/dE = N(r) as: 
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The l.h.s describes the spatial diffusion (K) and convec- 
tion (V c ), and the energy transport (first and second or- 
der terms). The r.h.s corresponds to the primary, sec- 
ondary and tertiary (only for antinuclci) source terms, 
and the sink (spallativc destruction) for the considered 
species. The most general form for the diffusion coeffi- 
cient is K(E, r). Energy gain and losses depend on (E, r) 
as well. The first order term f (E,r) corresponds to 
the sum of four contributions: ionization, Coulomb and 
adiabatic losses, and first order reacceleration. Ioniza- 
tion losses take place in the neutral interstellar medium 
(ISM), while the Coulomb ones in the completely ion- 
ized plasma, dominated by scattering off the thermal 
electrons lfj, [U 



The spatial dependence of these two 
terms is encoded in the distribution of the neutral and 
ionized gas. Adiabatic losses are due to the expanding 
wind and their spatial dependence is related to the gra- 
dient of V c - The last contribution to f Q (E,r) has the 
same origin as the second-order term s (E,r). Those 
come from the scattering of the charged particles off the 



turbulent magnetic fields of the Galaxy. As well as be- 
ing responsible for spatial diffusion, the Alfvenic waves 
also lead to energy drift and reacceleration. A minimal 
reacceleration scheme is well-motivated 12j and allows to 
calculate the /„ and s coefficients. Similar albeit more 
empirical forms have also been used 13j, [l4[ . In all these 
models, the strength of the reacceleration is mediated via 
the Alfvenic speed V a of the scatterers. 



A. The two-zone disk-halo model 

A full numerical treatment is generally required to 
solve the transpor t eq uation, as described, e.g., in Strong 

However, analytical (or semi- 
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and Moskalcnko 
analytical) solutions may be derived assuming a simpli- 
fied description of the spatial dependence of some param- 
eters in Eq. (fT]). The two-zone diffusion model 16l. Il7j|. 
based on the description of the Galaxy as a thin gaseous 
disk embedded in a thick diffusive halo, proved to be suc- 



cessful in reproducing the nuclear |14|. 18l . Il9j , antiproton 
[20I ] and radioactive isotopes data 19 1. It also allows to 



treat contributions from dark matter (or other exotic) 
sources located in the diffusive halo 21 , 2^, 23 1 , which is 
the aim of this Paper. We remind below the salient fea- 
tures of this model, which has been extensively detailed 

mRefs. mm nil. 

a. Geometry (L, R and h). The Galaxy is defined as 
a cylinder with a diffusive halo of half-height z — L and 
radius r = R. The halo thickness L is a free parameter 
of the model. The interstellar (IS) gas and the nuclei 
accelerators are contained in a thin-disk of half-height 
h <C L. The two parameters h and R are set to 100 pc 
and 20 kpc, respectively. 

b. Diffusion coefficient (Kq and 8). Diffusion arises 
because charged particles interact with the galactic mag- 
netic field inhomogeneities. The diffusion coefficient 
K(r, E) is related to the power spectrum of these inho- 
mogeneities, which is poorly known. Several analytical 
forms for K have been assumed in the literature. We 
consider here the standard rigidity (1Z = pc/Ze) depen- 
dent form K(E) = (3 Kq x 7Z s , where the normalization 
A'o is expressed in units of kpc 2 Myr -1 . The same diffu- 
sion coefficient is assumed throughout the Galaxy, i.e. in 
the disk and in the halo. Kq and 6 are free parameters 
of the model. 

c. Galactic wind and Alfvenic speed (V c and V a ). 
The convective wind is assumed to be of constant mag- 
nitude directed outwards perpendicular to the Galactic 
plane V c = V c e z . The reacceleration strength, mediated 
by the Alfven velocity V a , is confined to the thin disk. 
The first and second order terms f and s in Eq. (UJ) 
follow the formulation given in Ref. 14 1. 
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max 


0.46 


0.0765 


15 


5 


117.6 


39.98 


med 


0.70 


0.0112 


4 


12 


52.9 


25.68 


min 


0.85 


0.0016 


1 


13.5 


22.4 


39.02 



TABLE I: Transport parameters providing the maximal, me- 
dian and minimal primary antideuteron flux and compatible 
with B/C analysis (x| /c < 40) [3, [3. 



The free parameters of this propagation model — the 
diffusive halo size L, the normalization Kq and the slope 
5 of the diffusion coefficient, the value of the constant 
galactic wind V c and the level of reacceleration through 
the Alfvenic speed V a — were constrained from the study 
of the B/C ratio in Refs. 



14l |24| . In this Paper we take 



advantage of the parameters found in Ref. [l4[ and listed 
in Tabic [IJ Actually, when fitting to existing B/C data, 
a strong degeneracy of the transport parameters is ob- 
served, meaning that many sets of these five parameters 
are acceptable and lead to the same B/C ratio, but also 
to the same secondary (standard) antiproton flux [2p| . 
This degeneracy is broken for sources located in the dif- 
fusive halo, leading to large astrophysical uncertainties 



for the relevant fluxes [2l|, |26j, |27j ■ The same conclusions 
are obtained here for antideuterons. We will come back 
to this point in Sees. HVIandlVl 



B. The case of antideuterons 

Once the astrophysical framework for the transport of 
(anti)nuclci is set, the calculation of the antideuteron flux 
rests on the d specificities regarding the source term and 
its nuclear interactions [r.h.s. of Eq. (JTJ)]. The source 
term for antinuclei is usually cast in separate contribu- 
tions: 



1 (?,E)+Q Bec (f,E)+Q tei (r,E) 
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Among these three terms, only primary and secondary 
are true sources. They will be discussed in Sees. IIVI 
andlVl 



The tertiary term was emphasized in Ref. [28| to de- 
scribe the process corresponding to the non-annihilating 
interaction of a CR antinuclcus with an atom of the IS 
gas. In a medium of constant density n (v is the CR 
velocity): 
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The cross sections in Eq. <|3j> refer to inelastic non- 
annihilating processes and are detailed in the Appendix. 
The tertiary mechanism does not actually create new an- 
tideuterons. It merely states that the number of antin- 
uclei observed at energy E has to take into account the 
redistribution of those with energy E' > E (first term, 
positive contribution), minus the total number of d re- 
distributed to lower energies (second term, negative con- 
tribution). The tertiary contribution is treated as a cor- 
rective factor and dealt with iteratively: the equilibrium 
spectrum (r, E) is first calculated with Q^-j = 0, 
then q[H calculated with (r, E) in Eq. © to ob- 
tain N^-'{r,E), etc. For antideuterons, due to the small 
cross section — as shown in Ref. Q — only one iteration is 
necessary to converge to the solution (compared to a few 
iterations for antiprotons [13]). 

Hence, whether secondary or primary (or a mixture of 
both) sources are considered, three cross sections always 
enter the calculation: the differential non-annihilating 
inelastic cross section da non ~ ann /dE, the total non- 
annihilating inelastic cross section cr noI1 - ann anc j the to- 
tal annihilating inelastic cross section cr ann [appearing in 
the r.h.s. of Eq. (JTJ)]. Considering the ISM as a mix- 
ture of H and He, the full calculation requires six cross 
sections. Our calculations, based on the parameteriza- 
tions discussed at length in Duperray et al. 0, are dis- 
cussed in the Appendix where, in particular, we update 
d<r non ~ ann /dE. Note however, that even if some of these 
cross sections are slightly modified, the impact on the 
propagated spectra do not change the conclusions found 
in Ref. 0]. 

In the two-zone diffusion/convection/reaccclcration 
model, it is possible to extract semi-analytical solutions 
of Eq. ([T]), based on Bessel expansions of the transport 
equation. We do not wish to repeat the various steps of 
this derivation, nor to rewrite the complete form of the 
solutions, which have already been given in several pa- 
pers. The solution for the so-called antideuteron stan- 
dard source (of secondary origin in the galactic disk) 
and the numerical procedure to treat reacceleration is 
detailed in Ref. [2(j. The solution for an exotic source 
distributed in the whole diffusive halo of the Galaxy can 
be found in Ref. 2lj, l26j . Actually these two papers refer 
to p, but formally, the solutions apply to d as well. 



III. ANTIDEUTERON PRODUCTION 

The production of cosmic antideuterons is based on 
the fusion process of a p and n pair. One of the sim- 
plest but powerful treatment of the fusion of two or more 
nucleons is based on the so-called coalescence model 
which, despite its simplicity, is able to reproduce re- 
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markably well the available data on light nuclei and 
antinuclei production in different kinds of collisions. In 
the coalescence model, the momentum distribution of 
the (anti)deuteron is proportional to the product of the 
(anti)proton and (anti)neutron momentum distribution 
|29j, [30(. That function depends on the difference 
between (anti)nucleon momenta. It is strongly peaked 
around A^ ~ (compare the minimum energy to form a 

2.2 McV), so that 



4m p , with the binding energy 



kp 



k ~ kd 



(4) 



The d density in momentum space is thus written as the p 
density times the probability to find an n within a sphere 
of radius p$ around kp (sec, e.g. Ref. 31|): 
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The coalescence momentum po is a free parameter con- 
strained by data on hadronic production. Note that the 
coalescence model has been refined to account for heavy 
nuclei reactions [see, e.g. [32J, |33[, but as it is not relevant 



for this study, we will stick to the simple Eq. ([5]). 

The number dM^ of particles X produced in a single 
reaction R and which momenta are kx can be expressed 
as a function of the total available energy y/s, the inclu- 
sive (i.e. total inelastic or reaction cross section) and the 
differential cross section: 
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For instance, in our specific case X arc the antinucleons 
and antideuterons created in the pp, pHc and HeHe reac- 
tions between p-He CRs and H-He in the ISM. Assuming 
the usual equality between the unmeasured n and the 
measured p cross sections, and combining the two previ- 
ous expressions Eqs. ([5]) and ^ we get: 
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The hypothesis of factorization of the probabilities is 
fairly well established from experiments at high energies 
[see, e.g.,@[. For spallation reactions, however, the bulk 
of the antiproton production takes place for an energy 
y/s ~ 10 GeV, which turns out to be of the same order of 
magnitude as the antideuteron mass. Pure factorization 
should break in that case as a result of energy conser- 
vation. Two ansatz have been used in order to correct 
that effect for this regime: in Ref. P, Q it was assumed 
that, while the first antinucleon is produced with y/s, the 
center of mass energy available for the production of the 



second antinucleon is reduced by twice the energy car- 
ried away by the first antinucleon. Instead, in Duperray 
et al. Q the threshold production was phenomenologi- 
cally taken into account through an A + 2 phase space 
factor. The latter description seems more appropriate 
as, while preserving the correct asymptotic properties, it 
does not favor any mechanism for the pair production 
34 1. 

The coalescence momentum po is linked to the mea- 
sured coalescence factor B^ =2 (hereafter simply B 2 ): 
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The B 2 coefficient has been measured for proton-proton, 
proton-nucleus and heavy ion collisions (see a summary 
and references in Refs. [1, H^l). More recently, several 
other channels have also been measured at high energy: 
photo-production |35|. DIS production 36( and e + e~ 
production at the Z |3?1| and T(15) [11] resonances. The 
e + e~ channel is of particular interest for the DM annihi- 
lation reactions. 



A. Hadronic production 

For the hadronic processes, the coalescence momen- 
tum can directly be fitted to data. However, different 
assumptions regarding the set of data to retain can lead 
to different values of po. Note that many recent data are 
available for A + A systems, in addition to pp and pA 
reactions. However, the mechanisms at play in heavy ion 
collisions are not necessarily those of lighter systems (see, 
e.g., discussion in Sect. II. A of Ref. 0] and the results of 
Ref. [1^]), so that these reactions are discarded in the 
rest of our analysis. 

Chardonnct ct al. Q used pp data from Refs. 0, E], 
42j and pA collisions from Refs. 42j (see Fig. 1 in Ref. [Tj) 
to derive a coalescence momentum p = 58 MeV. Based 
on kincmatical relevance of the measured reactions, these 



authors disfavored the pp data from Ref. 4l| , but under- 
lined that a value po = 75 MeV, compatible with their 
whole set of data, would merely provide twice as more 
antideuterons. 

In Duperray et al. a larger set of data is used, 
including many pA reactions (see their Tab. I and ref- 
erences therein). The approach is more sound since a 
X 2 analysis on the momentum distribution of the frag- 
ments was performed, taking also into account the phase 
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space. This leads to an estimate of po = 79 McV. Note 
that at variance with the choice of Chardonnet et al. Q , 
Duperray et al. discarded the data from Ref. [42| be- 
cause they give a poor % 2 value compared to all other 
data. It is thus not surprising that these authors end up 
with a value close to p — 75 MeV quoted in Chardonnet 
et al. 0. In the present Paper, we take directly the cross 
sections derived in Duperray et al. @| using the value 
Po = 79 McV. 



B. Weak production 

At LEP energies, (anti)deuteron production occurs 
through e + e _ annihilations into qq pairs, a mechanism 
similar to the d production in DM annihilation reactions. 
Based on theoretical arguments, it has been argued 31 1 
that the antideuteron yields in e + e~ reactions should be 
smaller than in hadronic reactions. However, the ALEPH 



Collaboration [37f has found that this theoretical predic- 
tion (see Fig. 5 in ALEPH paper) underestimates their 
measured d inclusive cross section. They derive (see their 
Fig. 6) a value B 2 = 3.3 ± 0.4 ± 0.1 x 10~ 3 GeV 2 at the 
Z resonance, which translates into pa = 71.8 ± 3.6 MeV, 
very close to the po = 79 McV derived for the hadronic 
production. Hence, in the remaining of the Paper, the 
value of po = 79 MeV will be retained for both the pro- 
cesses of hadronic and electroweak origin. 



IV. SECONDARY ANTIDEUTERONS 

Secondary antideuterons are produced in the galactic 
disk from the collisions of cosmic protons and helium nu- 
clei over the ISM. We evaluate here the d propagated 
fluxes as well as the nuclear and propagation uncertain- 
ties, similarly to what was done for p in Ref. [2p| . 



A. Median flux 

The secondary d flux is the sum of the six con- 
tributions corresponding to p, He and p cosmic ray 
fluxes impinging on H and He IS gas (other reactions 
are negligible [8|). The p and He fluxes were fitted 
on BESS and AMS @, 0, 0] high energy data 
with a power law spectrum (see details in Ref. [2ol |') 
$(T) = N(T/GcV/n)~T. The best fit corresponds to 
N p = 13249 m~ 2 s" 1 sr- 1 (GeY/nT 1 and 7p = 2.72, 




and A Hc = 721 m" 2 s" 1 sr" 1 (GcV/n)- 1 and 7Hc = 2.74. 
The uncertainty on these two fluxes is small and leads to 
negligible uncertainties in the p and d spectrum. Contri- 
butions to the d flux from p + H and p + He reactions 
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FIG. 1: Contribution of all nuclear channels to the d sec- 
ondary flux. Dashed lines, from top to bottom refer to: p+H, 
p+He, He+H, He+He. Dotted lines, from top to bottom 
stand for: p+H, p+He. Solid line: sum of all the compo- 
nents. 



are evaluated using the p flux calculated in the same run. 
The production cross sections for these specific processes 
are those given in Ref. 0]. 

The different contributions to the total secondary an- 
tideuteron flux, calculated for the best fit propagation 
configuration (the "med" one in Table i.e. Kq = 
0.0112 kpc 2 Myr" 1 , L = 4 kpc, V c = 10.5 km s" 1 and 
V a = 52.1 km s _1 , are shown in Fig. [TJ As expected, 
the dominant production channel is the one from p-p col- 
lisions, followed by the one from cosmic protons on IS 
helium (p-He). As shown in Ref. t 8], the p+H channel 
is dominant at low energies, and negligible beyond a few 
GeV/n. The effect of energy losses, reacceleration and 
tertiaries add up to replenish the low energy tail. The 
maximum of the flux reaches the value of 2 ■ 10~ 7 parti- 
cles (m 2 s sr GeV/n)- 1 at 3-4 GcV/n. At 100 MeV/n it 
is decreased by an order of magnitude, thus preserving 
an interesting window for possible exotic contributions 
characterized by a flatter spectrum. 



B. Propagation uncertainties 

For a determination of the propagation uncertainties, 
we follow the same approach as in Ref. We calcu- 
late the secondary antideuteron flux for all the propa- 
gation parameter combinations providing an acceptable 
fit to stable nuclei 



14| | . The resulting envelope for the 



secondary antideuteron flux is presented in Fig. [5] The 



6 




1 10 

T s (GeV/n) 



FIG. 2: Dominant uncertainties on the interstellar secondary 
d flux. Solid lines: propagation uncertainty band. Dotted 
lines: nuclear uncertainty band. 



solid lines delimit the uncertainty band due to the degen- 
eracy of the propagation parameters: at energies below 
1-2 GeV/n, the uncertainty is 40-50 % around the aver- 
age flux, while at 10 GeV/n it decreases to ~ 15 %. This 
behavior is analogous to that obtained for p [2(3] and is 
easily understood. The degenerate transport parameters 
combine to give the same grammage in order to repro- 
duce the B/C ratio. Indeed, the grammage crossed by 
C to produce the secondary species B is also crossed by 
p and He to produce the secondary p and d. In short, 
a similar propagation history associated with a well con- 
strained B/C ratio explains the small uncertainty. With 
better measurement of B/C expected soon, e.g. from 
PAMELA or TRACER Q, this uncertainty will 
further decrease and could become negligible. 



as the value of po depends on the choice for da^. Hence, 
to be very conservative and to keep a simple approach, 
we have spanned all hadronic production cross sections 
in the range around their reference value. If we 

wished to translate this into an uncertainty on the coa- 
lescence momentum, this would lead to the effective range 
Po = 79^3 MeV. Finally, in order to estimate the maxi- 
mal flux with the most conservative attitude (the highest 
is the secondary flux the lowest is the chance to outline 
an exotic contribution), the non-annihilating cross sec- 
tion was doubled, as its value is probably only a lower 
limit (see Appendix). On the other hand, to evaluate the 
minimal flux, we switched off the p + H(He) — > d + X 
contributions, which intensity remains very uncertain. 

The dotted lines in Fig. [2] take into account the sum 
of all the possible uncertainties of nuclear source, as de- 
scribed above. At the lowest energies the flux is uncer- 
tain by almost one order of magnitude, at 100 GeV/n by 
a factor of 4. We have checked that the solar wind mildly 
decreases the IS flux at low energies but leaves the uncer- 
tainty magnitude unchanged. It is obvious from Fig. [2] 
that the uncertainties on nuclear and hadronic cross sec- 
tions (dashed lines) are more important than the ones 
coming from the propagation models (solid lines). We 
emphasize once more our conservative attitude in esti- 
mating the nuclear band. However, if no dedicated cam- 
paigns of measurements for these cross sections will be 
carried out in the future, these uncertainties are not likely 
to be significantly reduced. 



PRIMARY ANTIDEUTERONS 



The source term for primary d to be cast into Eq. ([T]) 



qf lm {r, z,E)=-qi (a ann v) 



dN d ( pDu(r,z) 



dE s 



(10) 



C. Nuclear uncertainty 

The possible nuclear uncertainty can arise from two 
different sources. The first one is directly related to 
the elementary production process dcr?. It was found in 
Ref. [i^l that this may be cast into a ±25% in the p prop- 
agated flux, so that it should be translated to a rough 
±50% in the d flux. Second, there is the uncertainty 
on the coalescence momentum po- Using an independent 
model (i.e. different from the coalescence scheme) for 
d production, Ref. Q found that, conservatively, the d 
background was certainly no more than twice the flux 
calculated with p — 79 MeV. 

To some extent, these two uncertainties arc correlated 



where (<r alm u)o is the thermal average of annihilation 
cross section times the WIMP velocity, dN^/dEg is the 
source spectrum, /9dm(?~, z) is the distribution of the DM 
in the Galaxy and m x is the WIMP mass. The quantity 
(canni')o depends on the particle physics model. If not 
differently stated, we fix its value to 2.3 • 10~ 26 cm 3 s _1 , 
which corresponds to a thermal CDM relic able to ex- 
plain the observed amount of cosmological dark matter 
49l . 5C, 51 1 . This will be our reference value for most of 
the analysis. The coefficient rj depends on the particle 
being or not self-conjugate: for instance, for a fermion it 
is 1/2 or 1/4 depending on whether the WIMP is a Ma- 
jorana or a Dirac particle. In the following we will adopt 
i] = 1/2. The quantity £ parameterizes the fact that the 
dark halo may not be totally made of the species under 
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scrutiny (e.g a neutralino or a sncutrino) when this can- 
didate possesses a relic abundance which does not allow 
it to be the dominant DM component (see e.g or 
}53l]). In this case p x = £pdm with £ < 1. For our ref- 
erence value for (a ann v)o clearly one has £=1. The DM 
candidate may then be identified with a neutralino 21 1 
or a sneutrino 5J| in various supersymmetric schemes, 



but for the purposes of our discussion it does not need to 
be specified. We in fact wish to maintain the discussion 
at the most general level: we just need to specify the final 
state particles produced in the DM annihilation process 
and the ensuing energy spectra. The final-state particles 
all belong to the Standard Model, and this allows us to 
perform our discussion on a totally general basis. We 
will at the end specify our candidate to be the neutralino 
and discuss experimental capabilities in the framework 
of some specific supersymmetric scheme. 

Below, we briefly recall the main steps for the calcu- 
lation of the source term, before focusing on the propa- 
gation of these antideuterons in the Galaxy (Sect. fV"C]) . 
which is one of the main novelty in this Paper. 

A. Antideuteron source spectrum 

The production of antideuterons from the pair- 
annihilation of dark matter particles in the halo of our 
Galaxy was proposed in [l| . The interest in this possible 
DM detection channel has been the physics case for the 
proposal of the GAPS experiment 3, i|,IH and it has also 
been considered in Rcfs. 

HE EH Iii3. 

As previously discussed (see Sect. IIIip the production 
of a d relies on the availability of a p — n pair in a single 
DM annihilation. In the case of a WIMP pair annihila- 
tion, the differential multiplicity for antiproton produc- 
tion may be expressed as 



dN p 
~d~Ep~ 



F.h 



dE p 



The annihilation into a quark or a gluon h is realized 
through the various final states F with branching ratios 
£T h • Quarks or gluons may in fact be directly produced 
when a WIMP pair annihilates or they may alternatively 
result from the intermediate production of Higgs bosons 
or gauge bosons. Each quark or gluon h then generates 
jets whose subsequent fragmentation and hadronization 
yield an antiproton energy spectrum dNp /dEp. 

As in Ref . [l[ , we assume that the probability to form 
an antiproton (or an antineutron) with momentum kp 
(k n ), is essentially isotropic: 



dNp 

dE- p {x + X 



Air kp Ep Tp(-\fs = 2m, Ep) 



Applying the factorization-coalescence scheme discussed 
above leads to the antideuteron differential multiplicity 

2 



dNg 
dEj 



4pp 
3 fcj 



1 dE e 1 ^ 



F.h 



E; 



(13) 

We assume, as discussed in Sect. IIII1 that the same value 
of the coalescence momentum po = 79 MeV holds as for 
hadronic reactions. 

The evaluation of the differential antiproton spectrum 
dN^/dEp follows the treatment of Ref. We refer 

to this paper for the details of the p spectra from all 
the annihilation channels. The resulting d source spectra 
from different final states are not directly shown here. 
Instead, we will provide examples of propagated d spectra 
for the various final states in Fig. [3] 



B. Dark matter halo profile 

The distribution of DM inside galaxies is a very de- 
bated issue (see e.g. Ref. [HI] for a brief highlight on re- 
cent results and relevant references). Different analyses 
of rotational curves observed for several types of galax- 
ies strongly favour a cored dark matter distribution, flat- 
tened towards the central regions (Ref. [i^ and references 
therein). On the other side, many collisionlcss cosmolog- 
ical N-body simulations in A-CDM models are now in 
good agreement among themselves 60j , but for the very 
central regions some resolution issues remain open. It 
has been recently stressed that asymptotic slopes may 



not be reached at all at small scales 61, 62, 63 



>es may 

Mm- 



However, it is not clear whether the central cusp is steep- 
ened or flattened when the baryonic distribution is taken 
into account (e.g. 66, 6t]|). For definiteness, we consider 



a generic dark matter distribution: 



( n ) Px = PcdmM = p e {^} j-j- 



(r Q /a) c 



a ^ (0--y)/ot 



+ (r/a) 1 



(12) 



(14) 

where r Q = 8 kpc is the distance of the Solar System from 
the galactic center. The spherical pseudo-isothermal and 
cored DM profile with (a,/3,7)=(2,2,0) will be the ref- 
erence in our calculations. The total local — Solar Sys- 
tem — CDM density has been set equal to p = 0.42 GeV 
cm -3 , the core radius to a=4 kpc. This value and the 
total local density may be varied in large intervals by 
maintaining good agreement with observations. The an- 
tideuteron flux is very sensitive to the local distribution 
of dark matter p Q , since it appears squared in the deter- 
mination of the flux, while it is less sensitive to the chosen 
dark matter distribution function (as was already under- 
lined in Ref. [2l| for p) . For completeness and for compar- 
ison, we also consider some of the profiles obtained from 
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FIG. 3: Antideuteron flux for a WIMP mass m x =100 GeV 
annihilating into different final states: solid (black) line refers 
to bb, dotted (red) to uu, short-dashed (blue) to WW, long- 
dashed (green) line to ZZ. The dot-dashed (magenta) refers 
to it and m x =200 GeV. The annihilation cross section (here 
and in the following figures) is fixed at the value: (cr ann 7j)o = 



FIG. 4: Antideuteron flux for WIMPs of m x =50 GeV. Dot- 
ted (black) lines refer to the interstellar flux, solid (red) lines 
stand for the top-of-atmosphere flux, modulated at solar min- 
imum. For each set of curves, the three lines refer to the max- 
imal, median and minimal propagation configurations defined 
in Table |U 



A-CDM simulations: i) a standard NFW profile having 
(a, (3, 7)=(1,3,1), with a = 21.746 kpc [68|, ii) the steeper 
DMS-1.2 (1,3,1.2) profile with a = 32.62 kpc 0, hi) 
and the modified NFW profile with an logarithmic slope 
(hereafter N04), with a = 26.4 kpc 6l| (similar to the 
Einasto profile [H). Scale radii for NFW and DMS-1.2 
profiles are taken from Ref. 70|], while the parameters 
for the N04 DM density distribution are the same as in 
Ref. 7l|. All these profiles are normalized to po = 0.42 



GeV cm , in order to isolate the effect of the local den- 
sity, which can be easily rescalcd in the flux evaluation. 
We do not include any boost factor due to halo substruc- 
tures. This conservative attitude is corroborated by the 
results of Lavalle et al. [58j , where it has been shown that 
the boost factor is typically close to unity: only for some 
extreme and unlikely configuration it can reach a factor 
of 10. 



C. Primary antideuteron flux and uncertainties 

In the present Section, we show our results for the 
propagated antideuteron flux from DM annihilation. We 
follow the prescriptions detailed in the previous Sections 
for the production and the propagation of antideutcrons. 
Top-of-atmosphere (TOA) fluxes are derived from the 
interstellar (IS) ones treating the effect of the solar mod- 
ulation with the force field approximation. If not differ- 



ently stated, TOA fluxes correspond to a solar minimum 
activity with modulation potential cj>=0.5 GV. For the 
reference propagation configurations, we refer the reader 
to Table [0 



1. Fluxes for various annihilation states 

Figure [3] displays the d flux from a WIMP of mass 
m x = 100 GeV. Each curve corresponds to different pure 
(i.e. with BR=1) annihilation final states: bb, uu, WW, 
ZZ and it (for which m x = 200 GeV). The aim of the 
figure is to show the effect on the observable flux of the 
different \-\ annihilation final states from which the d 
originate. The transport parameters are the "med" ones 
of Tab. [J (providing the best fit to B/C data) and the 
fluxes are not modulated (IS spectra) . In the low energy 
part of the spectrum - around and below 1 GeV/n - it 
turns out that these fluxes show quite similar shapes and 
comparable normalization when varying the final state. 
This energetic range is the one in which a primary flux 
might emerge from the secondary counterpart. In addi- 
tion, as we will also discuss at the end of our Paper, it is 
the window explorable by experiments in a near future. 
For these reasons, we will adopt the antideuteron yield 
from an annihilation into a pure bb final state as a simple 
but representative case for our discussions. 
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FIG. 5: Uncertainty due to propagation models on the an- 
tideuteron (black solid lines) and antiproton (red dotted lines) 
interstellar fluxes. The WIMP mass has been fixed at the 
value m x =50 GeV. For each set of curves, the three lines 
refer to the maximal, median and minimal propagation con- 
figurations defined in Table [T] 



2. Propagation uncertainties 

Fig.[5]shows the uncertainties on the primary d flux due 
to propagation parameters. The three curves (dashed 
line: IS fluxes; solid lines: TOA fluxes) correspond to 
the maximum, median and minimal set of propagation 
parameters as gathered in Table |H The band between 
the upper and lower curve estimates the uncertainty due 
to propagation. At the lowest energies of hundreds of 
MeV/n the total uncertainty reaches almost 2 orders of 
magnitude, while at energies above 1 GeV/n it is about 
a factor of 30. The figure refers to a WIMP mass of 50 
GeV but the results are insensitive to this parameter, as 
well as from the solar modulation. The magnitude of 
the propagation uncertainty is similar to the one affect- 
ing the primary antiprotons [2l| , as explained in Fig. 
This behavior is drastically different from that observed 
on the secondaries (see Fig. [2]). Indeed, their propagation 
history is very different. Whereas secondaries originate 
from standard sources in the thin disk of the Galaxy, ex- 
otic primaries are produced in all the diffusive halo of 
the Galaxy. As shown in Rcf. 



72j, these primary antin- 



uclei do not suffer large energy losses, reacceleration or 
tertiary redistribution as they rarely cross the thin disk. 
Most of them arrive at Earth — substantially unshiftcd 
in energy — from an effective diffusion cylinder of height 
L* = min(L, 2K/V C ) and radius of a few L* centered on 
the observer 73j. Hence, the parameters driving the un- 



FIG. 6: Effect of changing the DM halo density profile, for 
a m x =50 GeV WIMP and for the "max" (solid) and "med" 
(dotted) configurations of Table [I] The effect is shown as the 
relative change in the IS antideuteron flux as compared with 
the reference case of a cored isothermal profile. The lower 
(black) lines refer to the NFW profile 16811 . the median (blue) 
lines to a cuspy profile with 1.2 slope [69( and the upper (red) 
ones to the N04 profile of Ref. [gj]. 



certainty on the primaries are, at high energy, the range 
allowed for the halo size L (sec Tab|T]), whereas, at lower 
energy, this uncertainty is further increased by the effect 
of the galactic wind (2K/V C becomes smaller than L). 

When comparing in details the p and the d fluxes from 
Fig. [5j the following conclusions can be drawn. First, 
the antiproton fluxes are a factor of 10 4 higher than the 
antideuteron ones, as expected from the fusion process 
into d. Then, at high energies the difference between 
the two fluxes is to be ascribed to their source spectra, 
which for antidcutcrons is the square of the antiproton 
one. This effect, added to the different weight of destruc- 
tion cross sections, is visible also in the lower energy tail 
of the spectrum. The destruction of the antideuteron nu- 
clei on the ISM alters the flux by a factor of two, while 
the antiproton one is modified by a mere 20-25%. 



3. Dark Matter Halo profile uncertainty 

The effect of changing the DM distribution function 
Pdm{i~,z) on the d flux is demonstrated in Fig. [6j We 
only modify the shape of the density distribution (as 
discussed in Sect. IVB[) . while keeping frozen the local 
DM density to p o =0.42 GcV/cm 3 . The DM mass is 
m x =50 GeV, but as explained in the previous Section, 
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FIG. 7: Interstellar and Top-Of-Atmosphere (TOA) an- 
tideuteron fluxes. The dashed (blue) line shows the primary 
flux for m x =50 GeV and (cr arm u) = 2.3 • 1(T 26 cm 3 s"\ 
the (red) dotted line denotes the secondary component and 
the (black) solid line stands for the total (signal+background) 
flux. Propagation model is the median one in Table [I] 



the source term mostly factors out, so that these con- 
clusions hold for any neutralino mass. We plot the ratio 
(<^j— <#f f )/<#f f where </>™ f is the reference flux calculated 
with the cored isothermal profile and 4>d corresponds to 
the NFW, DMS-1.2 and N04 profiles (see Sect. EE]) . The 
two classes of curves correspond to the maximal (upper, 
solid) and median (lower, dotted) propagation parame- 
ters. The difference on the fluxes calculated with the 
minimal set of propagation parameters (not shown) is 
negligible. 

The increasing steepness of the profile in central re- 
gions of the Galaxy is responsible for an increasing of 
the d flux which is more relevant for higher diffusive ha- 
los. In the case L = 15 kpc, the d obtained with a 1.2 
cuspy profile [6^] is a factor of 2 higher than the cored 
one, while the NFW |60| halo gives fluxes 30-40% higher 
than the isothermal one depending on energy. The high- 
est flux is obtained with the log-slope NFW-like profile of 
Navarro et al. [6l|, which predicts, among the considered 
DM profiles, the highest DM density in a wide radial in- 
terval around the Solar System, although it is flatter than 
the DMS-1.2 and NFW profiles in the central kpc of the 
Galaxy. The flux obtained with the median parameters 
(L = 4 kpc) is less significantly modified by a change 
in the halo profile. Indeed, a charged particle produced 
around the galactic center can more easily reach the So- 
lar System when the magnetic diffusive halo is larger and 
when it is more energetic. 



FIG. 8: TOA fluxes for primary (solid lines) and secondary 
(dashed line) antideuterons for the median propagation pa- 
rameters. From top to bottom, the solid lines refer to WIMPs 



with mass m x =50, 100, 500 GeV. 



VI. POTENTIAL FOR DETECTION: RESULTS 
AND DISCUSSION 

We now turn to the determination of the total flux we 
can expect from the standard astrophysical source (see 
Sect. IIV[) added to a possible contribution from DM an- 
nihilation. In Fig. we show the IS and TOA (solar min- 
imum) secondary and primary d fluxes, and their sum. 
The primary flux is for a WIMP with mass m x = 50 
GeV, annihilation cross section (cr a nnw)o = 2.3 ■ 10~ 26 
cm 3 s _1 , an isothermal profile and = 0.42 GcV/cm 3 , 
as in the previous figures. The discrepancy between pri- 
mary and secondary flux for Xj < 2 GeV/n is striking. A 
signal from DM annihilation as the one in our example 
would definitely increase by a large amount the number 
of expected antideuterons in the lowest energy bins with 
respect to the purely secondary flux. At 100 MeV/n the 
expected d flux from a cosmologically dominant DM par- 
ticle of 50 GeV mass is two orders of magnitude larger 
than the secondary d flux calculated within the same 
propagation model. One has to remind that the primary 
flux scales as m 2 : this means that, in the low energy 
sector, the signal can overwhelm the background up to 
masses of the order of few hundreds of GeV. This figure 
demonstrates that the search for cosmic antideuterons 
is definitely one of the most powerful indirect detection 
means for the DM annihilation in the halo of our Galaxy. 
The discrimination power between signal and background 
can be as high as few orders of magnitude. A major limit 
to this kind of experimental inspection may reside in the 
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FIG. 9: Ratio of the primary to total (signal+background) 
TOA antideuteron flux. Solid (black) curve refers to a WIMP 
mass of m x =50 GeV and for the MED propagation parame- 
ters. Dotted (black) lines show the MAX (upper) and MIN 
(lower) cases. Dashed lines refer to the MED propagation 
parameters and different masses, which are (from top to bot- 
tom): m x = 10, 100, 500 GeV (red, blue, magenta respec- 
tively) . 



tiny level of the expected flux (about four orders of mag- 
nitude less abundant than antiprotons), which neverthe- 
less is foreseen to become experimentally accessible in the 
near future i, 0, 0, III • 

Figure [8] displays the TOA d flux for the median propa- 
gation parameters and at solar minimum. Together with 
the secondary flux, we plot the primary one for three dif- 
ferent WIMP masses: 50, 100, 500 GeV and for the same 
reference value of the annihilation cross section. As dis- 
cussed above, lighter WIMPs would provide a striking 
signal, and sensitivity is present for masses up to few 
hundreds of GeVs. 

The discrimination power between primary and sec- 
ondary d flux may also be deduced from Fig. [9] The 
ratio of the primary to total TOA d flux is plotted as a 
function of the kinetic energy per nucleon, for the three 
representative propagation models and different WIMP 
masses (the annihilation cross section is again fixed at 
the reference value). This ratio keeps higher than 0.7 for 
Tj < 1 GeV/n except for m x =500 GeV. For propaga- 
tion models with L > 4 kpc - which is a very reasonable 
expectation - this ratio is at least 0.9 for masses below 
100 GeV. Increasing the WIMP mass, we must descend 
to lower energies in order to maximize the primary-to- 



secondary ratio. However, for a m x =500 GeV WIMP 
we still have a 50-60% of DM contribution in the 0.1-0.5 



FIG. 10: TOA primary (red solid lines) and secondary (black 
dashed line) antideuteron fluxes, modulated at solar mini- 
mum. The signal is derived for a m x =50 GeV WIMP and 
for the three propagation models of Table [T] The secondary 
flux is shown for the median propagation model. The upper 
dashed horizontal line shows the current BESS upper limit 
on the search for cosmic antideuterons. The three horizon- 
tal solid (blue) lines are the estimated sensitivities for (from 
top to bottom): AMS-02 [zl, GAPS on a long (LDB) and 
ultra-long (ULDB) duration balloon flights [1,|1,E1- 



GeV/n range. Of course, the evaluation of the theoreti- 
cal uncertainties presented in this Paper must be kept in 
mind while confronting to real data. Fig. [9] clearly states 
that the antideuteron indirect DM detection technique is 
probably the most powerful one for low and intermediate 
WIMP-mass haloes. 

We finally discuss in Fig. [TU] a possible experimental 
short term scenario. The secondary d flux for the me- 
dian configuration of Table |l] is plotted alongside the pri- 
mary flux from m x =50 GeV, calculated for the max- 
imal, median and minimal propagation scenarios. The 
present BESS upper limit on the (negative) antideuteron 
search Q is at a level of 2-10" 4 (m 2 s sr GeV/n)" 1 . We 
also plot the estimated sensitivities of the gaseous an- 
tiparticle spectrometer GAPS on a long duration balloon 
flight (LDB) and an ultra-long duration balloon mission 
(ULDB) @, pUll, and of AMS-02 for three years of 
data taking [74|. The perspectives to explore a part of 
the region where DM annihilation arc mostly expected 
(i.e. the low-energy tail) are very promising. If one of 
these experiments will measure at least 1 antideuteron, 
it will be a clear signal of an exotic contribution to the 
cosmic antideuterons. Note that for AMS, a sensitivity 
at the level of the one at low energy should be obtained 
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FIG. 11: GAPS ULDB reach compared to predictions for 
neutralino dark matter in low-energy supersymmetric mod- 
els, shown in the plane effective annihilation cross section 
C 2 (°"annt , )o vs. neutralino mass m x . The solid, long-dashed 
and short-dashed lines show our estimate for the capability of 
GAPS ULDB of measuring 1, 10 and 100 events, respectively, 
for the median propagation model of Table[I] The scatter plot 
reports the quantity £ 2 {o- alm u)o calculated in a low-energy 
MSSM (for masses above the vertical [green] dashed line) and 
in non-universal gauging models which predict low-mass neu- 
tralinos 0, EH, Sill El, El & [Red] Crosses refer to 
cosmologically dominant neutralinos, while [blue] dots stand 
for subdominant neutralinos. Grey point are excluded by an- 
tiproton searches. 



beyond 2.3 GeV/n thanks to the RICH [7J, |8l|- This 
higher energy region would be complementary to a low 
energy detection. 

Prospects for antideuteron searches for specific dark 
matter candidates are shown in Figs. [11] [12] and [T3l 
where the expected capabilities of GAPS ULDB are con- 
fronted with theoretical predictions for neutralino dark 
matter in various supersymmetric schemes. The results 
are shown in terms of the effective annihilation cross sec- 
tion ^ 2 ((Tann'u)o as a function of the neutralino mass. We 
recall that DM annihilation signals are proportional to 
the square of the relic-particle density profile p x (which is 
in turn proportional to the total DM density through the 
rescaling factor £ defined in Sect.fV}, to the low-energy 
thermally-averaged annihilation cross section (cr ann u)o 
and inversely proportional to the square of the WIMP 
mass m x . It is therefore convenient to show our pre- 



FIG. 12: The same as in Fig. 1111 except that the supersym- 
metric predictions refer to a minimal SUGRA scheme. 



mate for GAPS ULDB to measure 1, 10 and 100 events, 
respectively. The number of events has been obtained 
by using an acceptance x time for the GAPS detector in 
its whole energy range of 1.33 x 10 7 m 2 sr s GcV/n 
for which no contribution to the antideuteron signal is 
expected from the secondary component. The quoted 
number of events therefore refers uniquely to the signal 
component. The transport parameters used in the de- 
termination of the antideuteron fluxes are those which 
provide the median flux and a minimum of activity for 
the Sun has been considered for the solar modulation. 
For definitcness, the event-lines have been calculated for 
the case of a pure bb DM annihilation final state. Fig. [3] 
shows that no major differences are present when differ- 
ent final states are considered: we therefore consider the 
lines of Fig. [TT] as a good estimate of the GAPS ULDB 
sensitivity, to be confronted with theoretical predictions 
for specific candidates. 

The scatter plot shows the distribution of points ob- 
tained in two supersymmetric models: for masses larger 
than about 50 GeV, the points refer to a low-energy Min- 
imal Supersymmetric Standard Model (MSSM), where 
all the mass parameters are defined at the electroweak 
scale 



21 



diction in the plane £ (cr an n^)o 



The solid, long- 



dashed and short-dashed lines in Fig. [TT] denote our csti- 



| . For masses lighter than 50 GeV, the scatter 
plots refers to the fight-mass neutralinos in the models 
of Refs. @, EH, EE E3, El, El HI where the gaug- 
ino universality condition is violated. The (red) crosses 
denote configurations for which the relic neutralino is 
a dominant dark matter component (i.e. its calculated 
relic abundance falls in the WMAP range 0, [ 



51] 
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FIG. 13: The same as in Fig. 1111 except that the astrophysical propagation parameters are those which predict minimal (left 
panel) and maximal (right panel) antideuteron fluxes. In the right panel, grey point are excluded by antiproton searches. 



for the CDM content of the Universe), while the (blue) 
points refer to subdominant neutralinos (for these points 
£ = O x /i 2 /(^cdm^ 2 )wmap applies). In the case of the 
standard MSSM, our prediction is that up to about a 20 
events could be detected, for masses between 50 and 100 
GeV. In the case of the low-mass neutralino models, a 
large number of events (up to a few hundreds of events 
for 10 GeV neutralinos) are expected for most of the rele- 
vant supersymmetric parameter space. In this situation, 
where the antideuteron fluxes are large, also antiproton 
production is expected to be sizeable. We therefore need 
to check that we do not produce an excess of exotic an- 
tiprotons which goes in conflict with current measure- 
ments on this antimatter channel. We employ here a 
2-cr upper bound on the admissible excess of antiprotons 
over the background in the lowest energy bin of BESS 



M185J and AMS [86( of 0.68 x lO^m-^r^s^GeV 



These configurations, which exceed the antiproton bound 
771 ] . are shown by the gray points. We see that, once 
the p bound is included, up to about 100 d events are 
expected for neutralinos of 15-20 GeV. A detector like 
GAPS ULDB will therefore have a strong capability to 
access this class of supersymmetric models, which are 
currently only partially probed by antiprotons searches 
771 ] and direct detection [87|, depending on the actual 
galactic parameters, like profile and local density. In the 
case of the standard MSSM, GAPS ULDB sensitivity al- 
lows to probe neutralino masses up to 300 GeV. 

The case of supersymmetric models endowed with 
Supergravity (SUGRA) boundary conditions [g^ are 
showed in Fig. [T^] for the same astrophysical models of 



Fig. [TTJ together with our GAPS ULDB expectations. 
In the case of minimal SUGRA we predict up to a few 
events, for some specific supersymmetric configurations, 
although the bulk of the parameter space requires larger 
sensitivities to be probed. 



The cases of the minimal and maximal set of transport 
parameters is shown in Fig. [T31 again for the MSSM and 
gaugino non-universal models. In the case of the "min" 
set of parameters, light neutralinos may produce up to 
70 events, while the in the standard MSSM at most a 
few events are predicted. For the "max" case we again 
are confronted with the fact that an excess of exotic an- 
tiprotons in conflict with current p measurements plays a 
role. Once this bound is properly taken into account, we 
obtain that 100-300 antideuterons are expected in the 
low/mid-mass range for neutralinos. Also in the stan- 
dard MSSM model, up to 40 events are foreseen. In this 
case, the sensitivity of GAPS ULDB extends up to neu- 
tralino masses of the order of 500 GcV. 



In conclusion, antideuteron searches offer very good 
opportunity of detecting a signal, for a variety of super- 
symmetric schemes. Especially low mass neutralinos in 
gaugino non- universal models [H, EO, 0, [z3, H H 
provide expectations for a large number of events in a 
detector like GAPS ULDB, with good opportunities to 
clearly disentangle a signal. 
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VII. CONCLUSIONS 

In this Paper we present a novel and updated calcu- 
lation of both the secondary and primary antideuteron 
fluxes, with special attention to the determination of the 
uncertainties of nuclear and astrophysical origin which 
affect the theoretical predictions of the d flux. The galac- 
tic environment is treated in a two-zone diffusion model, 
the same that successfully reproduces cosmic-ra y nu clear 
data, like e.g. the boron-to-carbon ratio [3, 13, and 
that is able to predict the observed antiproton flux 21 1. 
We therefore refine our calculation of Ref. [l[ , where the 
antideuteron signal as a promising tool for dark matter 
searches was originally proposed. 

We review the nuclear and astrophysical uncertainties 
and provide an up to date secondary (i.e. background) 
antideuteron flux. Propagation uncertainties for the sec- 
ondary component range from 40-50 % around the aver- 
age flux at energies below 1-2 GeV/n down to ~ 15 % at 
10 GeV/n. Nuclear uncertainties are largely dominant: 
a generous and conservative factor of 10 at very low en- 
ergies, which reduces to a factor of four at 100 GeV/n. 

The primary contribution has been calculated for 
generic WIMPs annihilating in the galactic halo, for 
the different production channels through which the an- 
tideuteron signal may be produced by DM annihilation. 
As for the antiproton 21 1 and positron signals [89||, we 
obtain that the transport-related processes induce the 
largest source of uncertainty in the antideuteron flux: it 
ranges from a factor of 10 upward and downward with 
respect to the median, best-fit, prediction, in the low 
energy range. 

Our theoretical predictions have then been confronted 
with the expected sensitivities of future detectors, specifi- 
cally GAPS US, S3] and AMS HlzlSi]- We have con - 
sidered neutralinos as dark matter candidates and dis- 
cussed three specific supersymmetric scenarios: a low- 
energy MSSM scheme, a gaugino non-universal super- 
symmetric model and a minimal SUGRA framework. For 
these classes of models we have analyzed the potential- 
ity of the GAPS detector in a ultra long duration bal- 
loon flight and found that this detector will have the 
capabilities to detect up to a few hundred events for 
low-mass neutralinos (in the mass range from 10 to 50 
GeV) in the gaugino non-universal models, and up to 
tens of events for 50-100 GeV neutralinos in standard 
low-energy MSSM. The sensitivity of GAPS ULDB on 
the neutralino mass extends up to 300-500 GeV in the 
low-energy MSSM. 

In conclusion, antideuterons offer an exciting target for 
indirect dark matter detection for low and intermediate 
WIMP masses and future experiments will have a unique 
opportunity to clearly identify a signal. 
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FIG. 14: Modeling of the total (elastic+inelastic) cross sec- 
tion cr^ f d . All curves are displayed as a function of the mo- 
mentum. The black dashed line corresponds to the Glauber 
approximation used in Duperray et al. Crosses refer to 
data for the charge conjugate reaction p+d The black 
solid line is a crude fit of the latter reaction when data are 
available (the high energy regime comes from Ref. [9l[). The 
two dotted curves illustrate a tentative estimate of a dH as 
2a pp (blue) and <j pp + a pn (red) , where we have modelled the 
latter cross sections by fits on the available data [ill ]. 
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APPENDIX A: CROSS SECTIONS 

We review here the cross sections (total, elastic, inelas- 
tic and non-annihilating) used in the calculation for pri- 
mary or secondary d. Their detailed and thorough study, 
including the production, has been presented in Duper- 
ray et al. 0. We do not repeat their arguments, but 
rather recall the main characteristics of these cross sec- 
tions, providing alternative formulations, minor correc- 
tions or further checks when possible. However, none of 
these "updated" cross sections lead to significant changes 
in the calculation of the d flux, compared to that calcu- 
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FIG. 15: Elastic cross sections as a function of the inci- 
dent particle momentum. Yellow crosses and blue circles are 
CERN data [9l|] for p+p (the dotted line is a fit to these data) 
and p + n reactions. The elastic cross section cr^ is modelled 
using 1o v J (solid black line). The dash-dotted grey line cor- 
responds to the Glauber cross section used in Duperray et al. 
[81. See text for comments. 



FIG. 17: Differential redistribution cross section d+p —> d+X 
as a function of the d kinetic energy per nucleon. Solid curves 
are for the limiting fragmentation hypothesis as used, e.g., in 
Tan & Ng [92| and dashed curves are for the Anderson et al. 
parameterization (9^|. The four sets of curves correspond to 
four incident d energy: 0.55 GeV/n (black), 1 GeV/n (blue), 
2 GeV/n (red) and 4 GeV/n (grey). 
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FIG. 16: From top to bottom: total, inelastic, elastic and 
inelastic non-annihilating cross sections for d+H as a function 
of the d kinetic energy per nucleon. 



lated with the cross sections presented in Ref. [8[. The 
main source of uncertainty remains the production cross 
section discussed in Sect. 



1. Total and elastic cross section 

In Duperray et al. Q , the Glauber approximation was 
used. No data exist for the total cross section of the pro- 
cess d + H, but there are measurements for the charge 
conjugate reaction d + p. Figure [14] shows various mod- 
eling of cfo t . Along with the charge conjugate measure- 
ments and the associated fit function, wc also plot the 



pn data. In this Paper, we assume <j. 



dp 



cr+Z. , so that 



the total cross section is given by the black solid line of 
Fig. [H Note that the combination 2<jPp and a™ + u pn 
are very close to a d P data, overshooting it by a mere 10%. 

Concerning the elastic cross section, the situation is 
worse since no data exist even for the charge conjugate 
reaction. However, as seen in Fig. [TJJ the total cross 
section a dp is well approximated by 2a pp or a pp + a pn , so 
this should be the also case for the elastic cross section. 
Since data for pn are scarce (Fig. [T5] blue circles) as 
compared to those for pp (yellow crosses), we choose to 
approximate <J p ^ d f=a 2crJ\ This is shown on Fig. [15] as a 
black solid line, which is a factor of two larger than the 
Glauber description (black dashed line) used in Ref. 0. 



Inelastic, total non-annihilating and differential 
redistribution cross sections 



The inelastic cross section is simply obtained from: 

finel = Otot — CTcl j (Al) 

and is itself further expressed as: 

(A2) 



_ non — ami , ^ann 

Cinel = 0- inel + C7;„ 



The non-annihilating (non-ann) part corresponds to an- 
tideutcrons that interact inelastically with H, but survive 
the collision, loosing a fraction of their initial energy. 

In Ref. [l[, the tertiary contribution was neglected, ar- 
guing that d (or symmetrically d) incident on a nucleon or 
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on a nuclear target should have a small non-annihilating 
cross section owing of the small d (d) nuclear binding en- 
ergy. As discussed in Ref. [§] , this intuitive argument can 
be invalidated both on empirical evidence and on formal 
grounds (see discussion and references therein) . We stick 
to the empirical approach by using the cross section given 
in Ref. Q . The total inelastic non-annihilating cross sec- 
tion c^tCnon-ann ( N ^R in Ref. @]) is obtained by sum- 
ming up the p + d — > (mr)pd cross sections which are 
experimentally available. As seen in Fig. [16] (black solid 
line), the cross section peaks at w 4 mb. As no attempt 
was made to evaluate the (expectedly small) contribu- 
tions of the channels not known experimentally, the over- 
all evaluation is a lower bound 0]. In the same Fig. [TBI 
for reference purposes, we show all the cross sections used 
in this Paper. 

The last issue is the energy redistribution of the non- 
annihilated d. In most of the papers dealing with antin- 
uclei [e.g. Q], [2(| H3i the energy distribution of the sur- 
viving particles is based on the limiting fragmentation 
hypothesis [see references in@]: 



da' 



dH^dX 



T dH 



dE k 



meLnon— ann 



(A3) 



However, experimentally, the pp — ► pX differential cross 
section was shown to be largely independent of the lon- 
gitudinal momentum pf of the produced particles in the 
center of mass and can be written in the laboratory frame 
as 931: 



d 2 a(pp — > pX) 
dpdQ 

„2 



p j(E — /3pcost 
2irp t E~ 



61 Op 2 



exp 



Pi 



0.166 



(A4) 



where 7 and (3 are the usual Lorentz factor and particle 
velocity, and pt the transverse momentum of the particle. 
Following Duperray et al. the same functional form 
is used for d and we define: 



do 



Anderson 



dE k 



2?r j2 



d 2 cr(pp — ► pX) 
dpdfl 



sinddO. 



(A5) 



Note that an incorrect normalization was applied in 
Ref. Q for the final cross section. The differential cross 
section for the tertiary term reads: 



da 



dH^rdX 



dE k 



dH 

incl,non— ann 



" ^^.Anderson " 




r E 'k 


dE k 




Jo 



' Jk da 



Anderson 



-dE'i 



(A6) 



Note, however, that this correction does not change sig- 
nificantly the tertiary spectrum and conclusions given in 
Ref. [f| . The inelastic scattering spectra in the laboratory 
for 0.55 (black), 1 (blue), 2 (red) and 4 GcV/n (grey) d 
are shown in Fig . IT7l For each of these energies, the An- 
derson et al. 93| scheme more efficiently redistributes d 
at very low energy (a few hundreds of MeV), as compared 
with the other one. 

Finally, we have not yet discussed the d+He cross sec- 
tions. Here, for the total, elastic, non-annihilating cross 
section, we simply use the naive geometrical factor A 2 / 3 
to account for interaction on He gas in the ISM: 



(A7) 



The same factor is also assumed for the differential ter- 
tiary cross section. 
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